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کاربرد زمین آمار در شناسایی محل‌های با بتانسیل بالای فرار آب در سد تنگاب * 
مقاله پژوهشی 
9 وهی( یا تا 1۳ ری مظفری؟ امین عصاری( 
چکیده سد تنکاب در فاصله ٩۰‏ کیلومتری جنوب شیراز و ۱۲ کیلومتری شمال غربی شهر فیرو زآباد بر روی رودنحانه فیرو زآباد احداث شده است. این سد در 
زون زاگرس و بر روی تاقدی سآهکی پودنو قرار دارد. محا سبات بیلان ذشان داده ااست که در این سد احتمال فرا رآب از محزن و تکیهگاه‌های وجود دارد. در 
این پژوهش تحلیل اکتشافیمکانی داده‌ها شامل بررسی نقشه‌های موقعیت مکانی, پارامترهای آماری, نمودار فراوانی و تغییرنما بر روی متفیرهای شانعص 
کیفیت سنگ, لوژان و متغیرسیمان‌حوری صورت گرفته است و جهت شناسایی زون‌های تراوا از روش‌های شبیه‌سازی زمینآماری استفاده‌شده است. متغیرهای 
مذکور در ترازهای ارتفاعی محتلف, تغییرات زیادی داشته است. پیو ستگی مکانی سیمانحوری بر مبنای تغییرنمای بدون جهت, ۱۰۰ متر بد ست آمده است. در 
حالیکه اين پیوستگی در متغیرهای شاحص کیفیت سنگ و لوژان بر اساس تغییرنمای جهت دار» به ترئیب بین ۲۰۰-۶۰۰ و 80-۱۰ متر متغیر است. درنهایت با 
توجه به نشه‌های شبیه سازی سه متغیر لوژان. شاحص کیفیت سنگ و سیمان‌خوری و تجزیه و تحلیل نشه‌های احتمال و عدم قطعیت متغیرهای مذکور در 
منطقه موردمطالعه, مشاهده شد در دو محل واقع در مرکز محزن و شمال غربی سد, مناطقی با نفوذپذیری بالا وجود دارد. لذا پتانسیل فرا رآب نیز در این مکان‌ها 
بیشتراست. بررس ی گسل‌های منطقه موردمطالعه نیز نشان داد تراکم گسل‌ها بیشتر در جناح راست سد بوده و اين می‌تواند سبب افزایش فوذپذیری در این منطقه 


باشد. 


واژه‌های کلیدی سیمان‌خوری. شاحص کیفیت سنگ. لوژان مدل‌سازی زمین‌آمار. 


مقدمه داده‌های مکانی» بر اساس پارامترهای آماری به‌دست‌آمده 


نالا شاگاه مها دور بش وا ری و 
درزه‌ها و رفتار هیدروژئومکانیکی. یکی از ملزومات اساسی 
مراحل اکتشافی به شمار می‌رود. به‌دلیل پیچیده بودن شرایط 
زمین‌شناسی. اندازه‌گیری داده‌های هیدروژئولوژیکی مثل قابلیت 
انتقال و هدایت هیدرولیکی ممکن است عدم قطعیت قابل توجهی 
را نشان دهند. بر این اساس» تخمین و شبیه‌سازی داده‌های 
هیدروژئولوژیکی با روش‌های مختلفی انجام می‌شود. روش‌های 
زمین آمار به‌طور گسترده‌ای در تخمین پارامترهای هیدرولیکی 
آبخوان و بررسی عدم قطعیت پارامترهای هیدروژئولوژیکی و 
تقو نکنیکین لاهن 

یکی از روش‌های مدل‌سازی هندسه لایه‌های تراواء استفاده 
از روش شبیه‌سازی زمین آماری است. این روش در مدل‌سازی 
سیستم شکستگی‌ها به‌منظور مطالعه جریانات سیالات در درون 
سنگ‌ها استفاده می‌شود [1]. در این روش پس از تجزیه‌وتحلیل 


* تاریخ دریافت مقاله ۱۸۰۰/۸/۱۶ و تاریخ پذیرش آن ۱۶۰۰/۱۰/۱۹ می‌باشد. 


(۱) دانشجوی کارشناسی ارشد آب شناسی, دانشگاه سیستان و بلوچستان 


(۲) نویسنده مسئول: دانشیار گروه زمین شناسی, دانشگاه سیستان و بلوچستان. 


(۳) استادیار دانشکده زمین شناسی, دانشگاه تهران. 


() دکتری آب‌شناسی, دفتر تحقیق و توسعه. معدن سنگ آهن گهرزمین. سیرجان. 


شبیه‌سازی متغیرها انجام می‌پذیرد [2]. مطالعات متعددی در 
سطح ایران و جهان با استفاده از زمین آمار و مدل‌سازی در 
خصوص نفوذپذیری. زمین شناسی. هیدرولوژی ‏ و 
هیدروژئولوژی انجام‌شده است که ازجمله می‌توان به [3]. [44 
[45 [۰۱ ۷1 [۰۸ ۸9 [410 [11] »1121 141.1۳[ ۱۵ 
[۰22 [۰۱23 [۰124 [۰125 [26] و [27] اشاره کرد. 

تغییرنما اساسی‌ترین ابزار در زمین آمار است که برای تشریح 
ارتباط مکانی یک متغیر به کار می‌رود. تغییرنماه کمیتی برداری 
که درجه همبستگی مکانی و شباهت بین نقاط اندازه‌گیری شده 
را برحسب مربع تفاضل مقدار دونقطه و با توجه به جهت و 
فاصله آن‌ها نشان می‌دهد معمولاً تغییرنما را با علامت (ظ)7 نشان 
نی دهتل راز وانطی (۱) ماه ی شوه [8ه]: 


0( ((20- (ظ + )یج < ()۷ 


261۲ 5616866.1050 ۵) 211622 اعصهطه ز :۵11 


که در آن ()7 تعداد جفت نمونه‌های به‌کاررفته در محاسبه 
که در فاضله 1 از یکذیگر قراز دارند؛ (262 مقدار مشاهده شد 
متغیر موردنظر و (201+0 مقدار مشاهده شد متغیر موردنظر که 
با فا ان( فرار دار انش ری با خه ی قرف انش ره 
۶ اگر در محاسبه تغییرنما به ازای هر لگ معین. بردار ظ 
نظیر آن لگ در فضا در هر جهتی بتواند قرار گیرد. در این صورت 
وتان با هزات ماع کون یت مت ناف ۲۰۱] 
این گونه تغییرنماها برای بیان تغییرنمای کلی منطقه مناسب 
می‌باشند و در مواردی که تعداد نقاط اندک باشد این‌گونه 
تغییرنماها اطلاعات کلی‌تری در مورد منطقه به دست می‌دهند 
[۳۰]. پس ازآنکه یک‌روند و ساختار مشخص برای تغییرنمای 
بدون جهت محاسبه گردید. در مرحله بعدی می‌توان تغییرنمای 
جهت‌دار و ناهمسانگردی را محاسبه نمود. تغییرنمای جهت‌دار 
را می‌توان این گونه تعریف نمود. اگر در محاسبه تغییرنما فقط از 
زوج نقاطی که بردار ظ آن‌ها (برداری که دونقطه را به هم وصل 
شک ) وی اشتاد معنن او من کرار داش بافله انفاده کوی 
در این صورت آن تغییرنما را تغییرنمای جهت‌دار می‌گویند [۲4]. 
معمولاً در مطالعات زمین‌شناسی با توجه به هدف موردمطالعه که 
بیشتر اوقات بررسی ویژگی‌های زمین‌شناسی و به‌ندرت باهدف 
بررسی ویژگی‌های آماری است. نمونه‌ها ب‌صورت منظم برداشت 
نمی‌شوند. درنتیجه استراتژی‌های نمونه‌برداری سبب شده است 
که در مناطقی مقادیر بالای یک متغیر متمرکزشده باشند. هرچند 
این مناطق ممکن است اطلاعات مفیدی به دست بدهند ولی 
نماینده سایر نقاط باقیمانده نیستند. برای به دست آوردن یک 
تخمین مناسب از میانگین داده‌های جامع می‌باید نمونه‌های 
موجود به‌گونه‌ای وزن دهی شوند که داده‌های کلاستره شده اثر 
بیش‌ازحدی را بر روی تخمین میانگین نگذارند .برای دکلاستره 
نمودن داده‌ها دو روش دکلاسترینگ چندضلعی و دکلاسترینگ 
سلولی وجود دارد. پیوستگی مکانی در بسیاری از متغیرهای 
موردبررسی در داده‌های علوم زمین وجود دارد. به‌طوری که بین 
داده برداشته‌شده در یک نقطه و داده‌های برداشته‌شده در نقاط 
دیگر که به فاصله معینی از یکدیگر قرار دارند. یکسری 
مشابهت‌ها و تفاوت‌هایی وجود دارد. دو داده برداشته‌شده در 
فاصله کم. عم مشابهت بیشتری نسبت به دو داده برداشته‌شده 
در فاصله بیشتر دارند. اگرچه بررسی تعدادی از نمودارهای 
پراکندگی برداری » توصیفی نسبتاً کامل را از پیوستگی مکانی یک 


نشریه مهندسی عمران فردوسی 


کاربرد زمی نآمار در شناسایی محل‌های با پتانسیل بالای فرا رآب در سد تنکاب 


متغیر ارائه می‌دهند. ولی این نمودارها حاوی اطلاعات کمی 
می‌باشند و نیاز به‌نوعی پارامتر خلاصه کننده دارند. سه تابع 
خلاصه کننده پیوستگی مکانی یک متغیر. ازجمله توابع ضریب 
همبستگی . توابع کوواریانسی و تغییر نماها . تمامی این توابع 
از یکسری پارامترهای آماری برای خلاصه نمودن نمودارهای 
پراکندگی برداری و ارائه آن» به صورتی که نشان‌دهنده تغییرات 
پیوستگی مکانی یک متغیر بافاصله و جهت باشد. بهره می‌برند. 
از بین آن‌ها تغییر نما مهم‌ترین ابزار برای بررسی پیوستگی مکانی 
بین داده‌هاست [۳۰]. 

در استان فارس. سد مخزنی تنگاب بخشی از آب آبیاری 
دشت فیروزآباد را تأمین می‌نماید. محاسبات بیلان توسط شرکت 
آب منطقه ای استان فارس و مطالعات قبلی نشان می‌دهد بخش 
زیادی از آب مخزن از طریق قسمت‌هایی از تکیه‌گاه‌های سد فرار 
کیت پمستیر که ناه شید اب ات که انطر 
زمین‌شناسی در چه مناطقی در اطراف این سد. پتانسیل فرار آب 
بالا استت: بتابراین.در انن بفوهتن.سخی شله است: فا با استفاده 
از بررسی پیوستگی مکانی داده‌های متغیرهای لوژان. شاخص 
کیفیت سنگ و سیمان‌خوری گمانه‌های اکتشافی درنهایت 
شبیه‌سازی‌های زمین‌آماری متغیرهای ذکرشده در منطقه 
موردمطالعه به این پرسش پاسخ داده شود. لازم به ذکر است 
پیش‌ازاین تحقیق. در محدوده‌ی مطالعاتی پژوهش‌های 
هیدروژئولوژی و مدل‌سازی نیز انجام‌شده است [۱۱]. در 
مطالعه‌ای احتمال فرار آب در سد تنگاب را با استفاده از 
ردیاب‌های رنگی بررسی کرده است. مشاهدات نشانگر وجود 
ارتباط هیدروژئولوژی بین ساختگاه سد و چشمه‌های محدوده 
تنگاب (روستای آتشکده) است. بر اساس رفتار رنگ نمودهای 
به‌دست‌آمده. این ارتباط از طریق یک سیستم کارستی تک 
مجرایی برقرار گردیده است. در پژوهشی دیگر [17]. در منطقه 
سد تنگاب. با استفاده از بررسی سری زمانی پارامترهای هیدرو 
شیمی چشمه‌ها؛ به شناسایی منبع چشمه‌ها و ارتباط آن‌ها با 
یکدیگر پرداعتند. بخشی از آبخوان توسط ۲۳ چشمه در سمت 
جنوب تخلیه می‌شود. یون‌های اصلی. هدایت الکتریکی, دما و 
17 این چشمه‌ها اندازه‌گیری شدند. این پارامترها نشان می‌دهد 
که منابع چهار چشمه, احتمالاً در قسمت‌های جنوبی قرار دارند 
و منابع دیگر چشمه‌ها هم در طرف‌های جنوبی و شمالی هستند. 


این بدان معنی است که همه این چشمه‌ها یک منبع مشترک و 
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رضا دهشیبی - رضا جهانشاهی - مرنضی مظفری - امین عصاری 


یک حوضه آبریز مشترک دارند [18]. از نتایج تست لوژان برای 
ارزیابی تخلخل موْثر و توسعه کارست در سد تنگاب استفاده 
کرده‌اند. در اين تحقیق رفتار آبخوان کارستی, نقش منطقه اپی 
کارستی نقش تغذیه در آبخوان‌های کارستی و اثر آن بر هیدرو 
گراف چشمه. مورد ارزیابی قرار گرفت و ارتباط کارستی شدن 
بر میزان نفوذپذیری (لوژان) بررسی‌شده است. در اين پژوهش 
مشخص شد که یکی از مژثرترین عوامل تأثیرگذار بر هیدرو 
گراف و تغییرات نوسان آب. طول مجرای اصلی کارستی است 
[ ۱۹ ]. نشت در جناح راست سد تنگاب را موردبررسی قرارداده 
اند. این تحقیق بر مبنای مدل عددی با استفاده از نرم‌افزار 2 
درزمینه جریان آب در خاک انجام‌شده است. نتیجه پیش‌بینی 
پژوهش نشان داده است عمده نشت از طریق مناطق کارستی و 
درزه و شکاف‌ها رخ‌داده است. 


مواد و روش‌ها 
ماطقه رب رونطا لیر 


( ۱۲ کیلومتری شمال غربی شهر فیروزآباده در مسیر فیروزآباد - 
کوار) بر روی رودخانه فیروزآباد واقع‌شده است (شکل -۱) اين 
سد در زون زاگرس چین خورده و بر روی تاقدیس پودنو قرار 
دار این تافقسن عمدها از سارت آهکن آسمازی (الیگومیوسه) 
تشکیل‌شده است و از بالا به‌وسیله سازند مارنی رازک و از زیر 
به‌وسیله سازند مارنی پابده گورپی (پاللوسن-الیگوسن) 
محصورشده است [۲۰]. سازند رازک در محدوده مخزن گسترش 
دارد و به دلیل اینکه از مارن و تبخیری‌ها تشکیل یافته‌اند از 
قابلیت نفوذ و آبگذری قابل توجهی برخوردار نیستند. لذا تنهاه 
محدوده سازند آسماری که دیواره دره ساختگاه سد را می‌سازد 
قابلیت آبگذری دارد. آهمک های سازند آسماری از استعداد 
کارست شدگی بالایی برخوردار بوده و وجود میان لایه‌های 
مارنی و مارن آهکی باعث توسعه کارست در اطراف این 
لایه‌ها شده است. این پدیده‌های کارستی در امتداد درزه‌ها و 
گسل‌های عمده که شیبی نزدیک به قائم دارند و توسط مواد 


ریزدانه پرشده‌اند توسعه بیشتری دارند [۲۱]. 


جت 91۵ 799 420 420 290 ۰ 296 126 ۰ ۵ 12046 
۰ 


شکل ۱ (الف) تصویر ماهواره‌ای از موقعیت سد. (ب) موقعیت سد تنگاب در استان فارس و 


(ج) نمایی از سد تنگاب. جهت نگاه به جنوب غرب 
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نشریهةُ مهندسی عمران فردوسی 


روش کار 
در اين پژوهش, مطالعات زمین آماری بر روی داده های لوژان, 
شاخص کیفیت سنگ و سیمان خوری چاههای اکتشافی در 
محدوده سد تنگاب انجام شده است که شامل موارد: تحلیل 
اکتشافی-مکانی داده‌ها. دکلاستره نمودن داده‌ها. بررسی 
هیستوگرام يا نمودارهای جعبه‌ای متغیرهای مختلف: بررسی 
تغییرات متغیرها درترازهای ارتفاعی مختلف. بررسی رابطه 
احتمالی بین متغیرهای مختلف با استفاده از نمودارهای 
پراکندگی. تهیه نقشه‌های تغییرنما و بررسی وجود و یا عدم 
وجود آنایزوتروپی احتمالی. مدل‌سازی زمین‌آمار (با استفاده از 
نرم افزار زمین آماری ایساتیس) شامل: ایجاد شبکه جهت انجام 
شبیه‌سازی گاوسیء انجام شبیه‌سازی گاوسی. تعیین حد 
آستانه‌های مختلف و بررسی تغییرات متغیرها در منطقه 
موردمطالعه. پردازش شبیه‌سازی‌های انجام‌شده و کمی سازی 

عدم قطعیت و در نهایت شناسایی زون‌های تراوا می باشد. 

به‌کارگیری اصل توالی تحت مدل چند گوسی توابع 
تصادفی را شبیه‌سازی گوسی متوالی می‌گویند. این الگوریتم برای 
شبیه‌سازی یک مشخصه واحد با استفاده از مقادیر واحد و سپس 
به دست آوردن اطلاعات ثانویه و شبیه‌سازی مشترک چندین 
مشخصه وابسته. معرفی‌شده است. این الگوریتم برای انجام 
شبیه‌سازی به داده‌های نرمال استاندارد نیاز دارد. برای شبیه‌سازی 
پیوستگی مشخصه 2 در آا گره در شبکه مشروط به مجموعه داده 

1( , ... ,1 < ۵ ,(2)1/0 است. مراحل شبیه‌سازی گوسی متوالی 

به‌صورت زير است: 

. اولین گام برر سی کردن صحت مدل چند گو سی است که 
نیاز ب‌رسم هیسستوگرام برای داده‌های اولیه و نرمال کردن 
داده‌ها و تبدیل 7-218 به ۲-۵26۵ با یک تابع توزیع تجمعی 
(]60) استاندارد نرمال است. 

۲ اگر مدل چند گوسی برای متغیر ا حفظ شد. در این حالت 
می‌توان شبیه سازی گو سی متوالی را به صورت مراحل زیر 
برای ۷-21۵ انجام داد: 
- ترسیم تغییرنما با استفاده از داده‌های نرمال شده. 
- انتخاب م سیر 5صادفی برای انجام شبیه سازی به‌طوری که 

هر گره شبکه فقط یک‌بار شبیه‌سازی می‌شود. 
- در هر گره "له پارامترهای میانگین و واریانس تابع توزیع 
تجمعی گوسی تعیین می‌گردند و سپس با انجام کریجینگ 


نشریه مهندسی عمران فردوسی 


کاربرد زمی نآمار در شناسایی محل‌های با پتانسیل بالای فرا رآب در سد تنکاب 


ساده با استفاده از مدل تغییرنمای داده‌های نرمال شدهه. 
تخمین داده‌ها بر روی شبکه صورت مین کیر: 
- بیرون کشیدن مقدار شبیه سازی شده از تابع توزیع تجمعی 
و اضافه کردن آن به مجموعه داده‌ها. 
- رفتن به گره بعدی در طول مسیر تصادفی و تکرار دو 
مره ی 
تکرار کلیه مراحل بالا تا اينکه آا گره شبیه سازی شوند و 
کل شبکه دارای عدد شود [28]. 


- شبیه‌سازی گوسی متوالی یک‌راه سریع و ساده برای انجام 
شبیه‌سازی زمین‌آماری است؛ زیرا برای مدل‌سازی گوسی 
در هر مکان نیاز به حل کردن تنها یک سیستم کریجینگ در 
آن مکان است. اما درعین‌حال توانایی نشان دادن پیوستگی 
مکانی برای طبقات مختلف مقادیر را ندارد و همچنین 


نتایج و بحث 
تحلیل اکتشافی - مکانی 

بررسی آمار مقدماتی بر روی ۲۳۷۰ داده شاحص کیفیت سنگ. 
۷ داده لوژان و ۱۱۸۶ داده سیمان خوری گمانه‌های اکتشافی 
حفرشده در تکیه‌گاه‌های سد تنگاب انجام‌شده است؛ میانگین؛ 
واریانس و انحراف معیار متغیر شاخص کیفیت سنگ به ترتیب 
برابر با ۸۷/۳۸ ۳۰۶/۱۳ و ۸۰۸ ۱۹ است (جدول ۱). نمودار 
هیستوگرام فراوانی آن (شکل ۲) نشان‌دهنده چولگی منفی و عدم 
نرمال بودن داده‌ها است. ازآنجایی‌که دامنه تغییرات آن دارای 
مقادیر بالاست. این موضوع نشان‌دهنده وجود مقادیر بالا در کنار 
مقادیر پایین شاخص کیفیت سنگ در ترازهای مختلف ارتفاعی 
است. همان‌طور که اشاره شد داده‌ها از توزیع نرمال پیروی 
نمی‌کنند و جهت انجام شبیه‌سازی يا تخمین و محاسبات 
زمین آماری نیاز به نرمال کردن داده‌ها است. 

شکل (۲سب) نشان می‌دهد که میانگین متغیر شاخص 
کیفیت سنگ درترازهای مختلف (بین تراز های ۱۲۹۸ تا ۱6۵۰ 
متر از سطح دریا) داری نوسان است و ازآنجایی که دامنه تغییرات 
و پراکندگی داده‌ها زیاد است؛ لذا میانگین نمی‌تواند معیار خوبی 


جهت مقایسه این متغیر درترازهای مختلف باشد. 


سال سی و چهارم» شمارة چهار» ۱۶۰۰ 


رضا دهشیبی - رضا جهانشاهی - مرنضی مطفری - امین عصاری 


جدول ۱ آمار مقدماتی شاخص کیفیت سنگ. لوژان و سیمانخوری بدون تخصیص وزن 


۱۵۱۵ 4 


متغیر تعداد | کمینه | بیشینه | میانگین | انحراف معیار 
شاخحص کیفیت سنگ | ۲۳۷۰ ۰ ۱۰۰ ۸۷/۳۸ ۱۹/۰۸ 
لوژان 52 ۰ ۳۱۱/۳ 1/۱6۵ ۱۸/۰۹ 
سیمان‌خوری ۱۱۸۶ ۹ ۳۳۸/۲۰ ۱۸۳/۸ 
لف 


واریانس | چولگی | کشیدگی | دامنه تغییرات | ضریب تغییرات 
۱۳/۳ ۹/- 2/۵ ۱۰۰ ۰/۲ 
۳/۳۳ 5۹۹ 2۷۹ ۳۳۲ ۰/6۸ 
۸ ۳( ۱۳/۳۹۰ ۱۱۲۱۱۸۱ ۱/۹ 1۹۹ 


شکل ۲ هیستوگرام مقادیر شاخص کیفیت سنگ و تغییرات میانگین آن در ترازهای ارتفاعی مختلف 


در نمودار جعبه‌ای (شکل ۲-الف) تغییرات متغیر شاخحص 
بر کفت سنگ در ترازهای حز مختلف ارتفاعی نان می‌دهد مقدار 
میانگین شاخحص کیفیت سنگ در بیشتر ترازها از مقدار میائه کمتر 
ست؛ که این موضوع نشان‌دهنده چولگی منفی در اکثر ترازها 
بوده و متغیر به سمت مقادیر کمتر میل کرده است. در تراز 
رتفاعی ۷۰-۱۲۹۸/۰/ ۱۳۱۳ متر از سطح دریا مقدار میانگین 


بیشتر از مقدار میانه بوده که بیانگر چولگی مثبت در این تراز بوده 


ست. نمودار میله‌ای دامنه تغییرات در مقابل ترازهای ارتفاعی 
مختلف (شکل ۲-ب) نشان می‌دهد در تراز ارتفاعی ۱۳۷۶/۱۰- 
۰ د«متر از سطح دریا) کمترین دامنه تغییرات وجود دارد 
و اختلاف بین بیشترین و کمترین مقدار متغیر برابر ۱۰ است. در 
رتبه دوم تراز ارتفاعی /٩۰‏ ۱۲۹۸ تا ۱۳۱۳/۷۰ متر است که 
احتلاف مقادیر بيشینه و کمینه آن برابر ۶۰ است. وجود مقادیر 
بالا در دامنه تغیبرات؛ در بعضی از ترازهای ارتفاعی» پراکندگی 
شدید داده‌ها را در ترازهای مختلف نشان می‌دهد و به این 
معناست که در هر تراز مقادیر بالای متغیر می‌تواند در کنار مقادیر 
پایین قرار داشته باشد. لذا معمولا جهت مقایسه و بررسی تغییرات 
متغیر در ترازهای مختلف مقدار ضریب تغییرات برای هر تراز 


محاسبه می‌شود. هرچه مقدار ضریب تغییرات بیشتر باشد 


سال سی و چهارم. شمارة چهار» ۱۶۰۰ 


نشان‌دهنده تغییرات شدیدتر متغیر در آن تراز ارتفاعی است. با 
توجه آمار مقدماتی متغیر شاخص کیفیت سنگ در ترازهای 
مختلف بدون تخصیص وزن. مشاهده‌شده است که بیشترین 
مقدار انحراف معیار در تراز ارتفاعی ۱۳۲۸/۸۰-۱۳۱۳/۷۰ 
(مترازسطح دریا) و برابر با ۲۷/۵ است. با توجه به محاسبه 
ضریب تغییرات متغیر در ترازهای ارتفاعی. بیشترین تغییرات در 
تراز ارتفاعی ۱۳۲۸/۸۰- ۱۳۱۳/۷۰ (متر) قرار دارد و ضریب 
تغییرات آن برابر ۰/۳٩‏ است و کمترین مقدار ضریب تغییرات 
مربوط به تراز ۱۳۷۶/۱۰ تا ۱۳۸۹/۲۰ (مترازسطح دریا) است که 
مقدار آن ۵ است (شکل ۳-ج). همان‌طور که پیش‌ازاین اشاره 
شد با توجه به هیستوگرام متغیر شاخص کیفیت سنگ از توزیع 
نرمال پیروی نمی‌کند و بدین منظور, داده‌ها با روش تبدیل نرمال 
کاوسی به تزمال: کلیل‌شده اسیق. 

نتایجم آمار مقدماتی متغفیر لوژان بدون تخصیص وزن 
(دکلاسترینگ ) در جدول (۱) نشان می‌دهد که مقادیر ميانگین و 
چولگی به ترتیب برابر با 7/۱۵ و ۵/1٩‏ است. با بررسی میزان 
دامنه تغییرات برابر با ۲۱۱/۲ و شاخحص‌های پراکندگی مشخص 
است که مقادیر لوژان در ترازهای مختلف در نوسان بوده است 


و بیانگر تغییرات این متغیر در ترازهای مختلف است. ازآنجایی که 


نشریهة مهندسی عمران فردوسی 


دامنه تغییرات داده‌ها زیاد بوده و احتلاف میان بزرگ‌ترین و 
کوچک‌ترین داده لوژان در هر تراز مقدار بالایی را نشان می‌دهد؛ 
میانگین نمی‌تواند معیار خوبی جهت مقایسه لوژان در ترازهای 
متتاف باشد, کش رنه لزان در فران ارفاض تاف.فن شیر داز 
جعبه‌ای شکل (۶).نشان می‌دهد در بیشتر ترازها مقدار میانگین 
بزرگ‌تر از مقدار میانه است؛ که نشان‌دهنده چولگی مثبت در 
توزیع این متغیر است. مقدار دامنه تغییرات لوژان در برابر 
ترازهای ارتفاعی مختلف در نمودار میله‌ای شکل (۶-ب) نشان 
می‌دهد کمترین دامنه تغییرات لوژان در تراز ارتفاعی ۳6- 
۰ متر از سطح دریاء و بیشترین ضریب تغییرات در ترازهای 
ارتفاعی ۱۳۰۵۰-۱۲۰ و ۱۳۵۰ تا ۱۳۹۵ (متر از سطح دریا ) 
قرار دارد (شکل 4-ج). همچنین نتایج آمار مقدماتی میزان سیمان 
خوری در گمانه‌های اطراف سد (جدول ۱) . نشان می‌دهد توزیع 
این پارامتر به‌صورت نرمال نبوده و چولگی مثبت. نشان‌دهنده 
عدم تفارن داده‌ها است. 
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شکل ۳ نمودار جعبه‌ای (الف). دامنه تغییرات (ب) و ضریب تغییرات (ج 
شاحص کیفیت سنگ درترازهای ارتفاعی مختلف 


نشریهُ مهندسی عمران فردوسی 


کاربرد زمي نآمار در شناسایی محل‌های با پانسیل بالای فرا رآب در سد تنگاب 


همان‌طور که قبلاً بیان شد. عدم توزیع یکنواخت مکانی 
داده‌ها. به‌شدت پارامترهای آماری را تحت تأثیر قرار می‌دهد. به 
همین منظور لازم است قبل از تحلیل آماری داده‌هاء وزن مناسب 
برای هر داده تعیین گردد. در اين مطالعه از روش دکلاسترینگ 
سلولی در کلاستره نمودن متغیرهای شاحص کیفیت سنگ و 
لوژان استفاده‌شده است. مقدار میانگین متغیر شاخص کیفیت 
سنگ و لوژان قبل از دکلاسترینگ به ترتیب ۸۷/۳۸ و 1/۱۵ بوده 
است که پس از اعمال وزن دهی به ۹۱/۸۷ و ۵/٩۱‏ تغییر يافته 
ین اعمال وزن دکلاسترینگ بر روی‌داده‌ها تأثیر چندان 


شم 


زیادی بر مقدار میانگین تل‌استت 
یکنواخت داده‌ها است. 
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شکل ۶ نمودارهای جعبه‌ای متغیر لوژان. (الف) نمودارمیله‌ای دامنه تغییرات 


چاه 
4 
4 


ضریب تقییرات 


(ب) نمودارمیله‌ای ضریب تغییرات درترازهای مختلف(ج) 


سال سی و چهارم. شمارة چهار: ۱۶۰۰ 


رضا دهشیبی - رضا جهانشاهی - مرنضی مظفری - امین عصاری 


مدل‌سازی تغییرنما 
جهت انجام شبیه‌سازی و آنالیزهای زمین‌آماری نیاز به نرمال 
کردن داده‌ها و ترسیم تغییرنمای مقادیر نرمال شده و برازش 
مدل‌های مختلف بر تغییرنما تجربی مدل است. ازاین‌جهت در 
این مرحله پیوستگی مکانی بر اساس تغییرنمای بدون جهت و 
جهت‌دار مقادیر نرمال شده متغیر شاخحص کیفیت سنگ. لوژان و 
سیمان‌خوری موردبررسی قرارگرفته است. در شکل (۵-الف)؛ 
تغییرنمای بدون جهت متغیر شاخص کیفیت سنگ مشاهده 
می‌شود. دو مدل اثر قطعه‌ای و توانی به تغییرنمای تجربی برازش 
تا اس اه مدل گرا ۲۱/۹۸ و مل ال قطعه‌ای با مقدار 

۹ بالاترین سهم از واریائس کل را دارا است. تغییرنمای مقادیر 

نرمال شاخص کیفیت سنگ در سه جهت مختلف محاسبه و 
مدل‌سازی شده است که به شرح زیر است: 

۱. تغییرنمای مقادیر نرمال برای شیب و آزیموت صفر در شکل 
(۵-ب): با توجه به این شکل. مدل توانی به تغییرنمای تجربی 
برازش شده است. دامنه مدل توانی و دامنه مدل ۷/۳۰ متر 
است و مدل اثر قطعه‌ای با مقدار ۰/۷ بالاترین سهم از واریانس 
کل را شامل می‌شود. 
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شعل (0-ج): دو مدل اثر قطعه‌ای و توانی به تغییرذمای 
تجربی برازش شده است. سهم اثر قطعه‌ای از واریانس کل 
برابر ۵ سهم مدل توانی ۵ است و دامنه مدل توانی 
برابر ۶۲۷/۵۲ متر است. 

۳ تغییرنمای مقادیر نرمال برای شیب ٩۰‏ در شکل (۵-د): دو 
مدل اثر قطعه‌ای و توانی به تغییرنمای تجربی برازش شده 
است و اثر قطعه‌ای با مقدار ۰/٩۱‏ بالاترین سهم از واریانس 
کل را شامل می‌شود. مجموع دامنه مدل برازش يافته به این 
تغییرنمای تجربی برابر ۱۷/۸۷متر است. 

اشکال (1- الف و ب) تغییرنمای سه‌بعدی به ترتیب در 
صفحات 2۷۲ و 22 را ذ شان می‌دهد. از شکل مذکور می‌توان 
دریافت پیو ستگی مکانی در امتداد محور شرقی-غربی بیشتر از 

بوستگی سکانی در امقداد موی 2 رابهای کمانه‌ها استه: 

شکل زر تغییرنمای سه‌بعدی در صفحه ۲2 را نشان 

می‌دهد. با توجه به شکل مذکور پیوستگی مکانی در جهت ۷ 

بیشتر از جهت 2 است. همچنین سه مدل تغییرنمای جهت‌دار 

مقادیر نرمال شده متغیر شاخص کیفیت سنگ نیز با یکدیگر 

تلفیق و به‌صورت یک مدل واحد ارائه‌شده است (شکل ۱-د). 
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شکل ۵ (الف) مدل برازش یافته تغییرنمای شاخص کیفیت سنگ بدون جهت. (ب) تغییرنما در جهت آزیموت صفر و شیب صفر (ج) تغیبرنما برازش يافته در 


جهت آزیموت ٩۰‏ و شیب صفر (د) تغییرنما برازش شده در جهت شیب نود 


سال سی و چهارم. شمارة چهار» ۱۶۰۰ 
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- ق۱۱۵8 


اهوم ددع۳ 


شکل. 1 «الف) تغییرنمای سه‌بعدی شاخص کیفیت سنگ در صفحه ۰26۷ (ب) صفحه ۰262 (ج) صفحه ۷2 و (د) تغییرنمای تجربی به همراه مدل واحدبرازش 


یافته درامتدادمحورهای سه بعدی ناهمسانگردی 


تغییرنمای تجربی جهت‌دار و بدون جهت مقادیر نرمال 
متغیر لوژان به همراه مدل کروی توانی و کوبیک بر مدل تجربی 
برازش يافته در شکل (۷) اراثه‌شده است. همان‌طور که در شکل 
مشاهده می شود دامنه مدل توانی ۱۳۹/۸۰ متر و مدل اثر قطعه‌ای 
با مقدار ۰/۵۶ بالاترین سهم از واریانس کل را دارا است. همچنین 
نتایج محاسبه مدل‌سازی و تغییرنمای لوژان در سه جهت مختلف 
عبارت است از: 
تغییرنمای شیب و آزیموت صفر در شکل (۷-الف): دو مدل 
کرو آفر فطعاای به سای خحریی برازشن شاه یک داد 
تقل کپ ۲۲۷۱۱۲۵ یوار قطعهی مقیان غک 
بالاترین سهم از واریانس کل را شامل می‌شود. 
۵ تغییرنمای مقادیر نرمال برای شیب صفر و آزیموت ٩۰‏ در 
شکل (۷-ج): دو مدل اثر قطعه‌ای وکوبیک بر تغییرنمای تجربی 
براقشن شده اسست هم ار فقای او وارباتی کل یرای ۳ 
سهم مدل کوبیک ۰/۳۳است و دامنه مدل توانی برابر با 2*۳/۲۹ 


نشريهة مهندسی عمران فردوسی 


4و اف با 
1 تغییرنمای مقادیر نرمال در شیب ٩۰‏ در شکل (۷-د): دو مدل 
توانی و اثر قطعه‌ای به تغییرنمای تجربی برازش شده است. اثر 
قطعه‌ای با مقدار ۰/۵1 بالاترین سهم از واریانس کل را شامل 
می‌شود. مجموع دامنه مدل برازش یافته به اين تغییرنما برابر 
هي انس ربا اف مه مین کي ای اسی‌ت‌دان مقادیر 
نرمال لوژان. یک مدل واحد ارائه‌شده است. 

بعد از نرمال‌سازی داده‌های سیمان‌خوری. نمودار تغییرنمای 
تجربی مقادیر نرمال سیمان‌خوری, به همراه برازش مدل‌های 
کروی و توانی و کوبیک تهیه‌شده است (شکل ۸. بررسی 
تغییرنما در جهات مختلف نشان داد هیچ روندی در پیوستگی 
مکانی سیمان‌خوری مشاهده نمی‌شود و به همین دلیل تغییرنمای 
بدون جهت محاسبه گردیده است. با توجه به این تغییرنما 


پیوستگی مکانی این متغیر حدود ۰ متر بوده است. 


سال سی و چهارم. شمارة چهار: ۱۶۰۰ 
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شکل ۷ (الف) مدل برازش يافته تغییرنمای لوژان بدون جهت. (ب) تغییرنما در جهت آزیموت صفر و شیب صفر (ج) تغییرنما برازش یافته در جهت آزیموت 


۰ و شیب صفر و (د) تغییرنما برازش شده در جهت شیب نود 
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شکل ۸ تغییرنما برازش داده‌شده بر مقادیر نرمال سیمان خوری (مدل کروی» کوبیک و توانی) 


سال سی و چهارم. شماره چهار» 5 ۳۹ 
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01 
م‌ 


صصصاه- ۱ هطم 


: صوطوها و۷ 


نشريهة مهندسی عمران فردوسی 


۳۶ 


شبیه‌سازی زمین آماری 

با کمک تغییرنماهاء شبیه‌سازی گاوسی متوالی متغیرهای شاخحص 
کیفیت سنگ. لوژان و سیمان‌خوری. مدل‌سازی گردید. جهت 
نمايش این شبیه‌سازی به‌عنوان مثال شکل (۹-الف)» نمایی 
سه‌بعدی از بلوک شاخص کیفیت سنگ می‌باشد و شکل (۱۰- 
ب). حد آستانه در ترازهای مختلف آن مشخص‌شده است 
به‌گونه‌ای که فقط مقادیر بالا یا پایین از یک حد آستانه مشخص 
را در ترازهای مختلف نشان می‌دهد که در این شکل برای 
شاخص کیفیت سنگ مقدار کوچک‌تر از ٩۰‏ درصد درنظرگرفته 
شده است. 

بررسی مقادیر متغیرها در سه بعد و در تمامی خروجی‌های 
شبیه‌سازی کاری دشوار است و نیازمند به احاطه کامل بر تمامی 
نتایج شبیه‌سازی است؛ لذا لازم است تا شبیه‌سازها در هر سلول 
و یا هر منطقه معین به‌صورت خلاصه‌شده به نحوی تعیین شود. 
بدین سبب نتایج صد شبیه‌سازی هرکدام از سه متغیر شاخص 
کیفیت سنگ. لوژان و سیمان‌خوری برهم منطبق شدند و مقادیر 


الف 


7 


ارو و 25255 
(ها2 


کاربرد زمی نآمار در شناسایی محل‌های با پتانسیل بالای فرا رآب در سد تنگاب 


میانگین و احتمال کوچک يا بزرگ‌تر از مقادیر خاص برای هر 
سلول محاسبه و در نظر گرفته‌شده است. شبکه شبیه‌سازی‌شده 
متغیرها در راستای گمانه‌ها (جهت 2 شامل ۲۵ سلول با 
اندازه‌ای برابر با ۱۰ متر است که درواقع ۲۵۰ متر شبیه‌سازی در 
راستای گمانه‌ها انجام‌شده است. در هر تراز ارتفاعی» کمترین 
شاخص کیفیت سنگ و بیشترین لوژان و سیمان‌خوری در مناطق 
معینی متمرکزشده است که در این مناطق میزان خردایش و 
نفوذپذیری بیشتر و احتمال دارای تراوایی بیشتر خواهد بود. لذا 
می‌توان مکان‌های مناسب را جهت تزریق سیمان در مرحله 
تکمیلی جهت جلوگیری از فرار آب را پيشنهاد دهد. هرچند در 
این نقشه‌ها مقادیر احتمال متفاوت است؛ اما روند کلی تغییرات 
در تمامی این نقشه‌ها یکسان خواهد بود. به‌عنوان‌مثال میانگین 
تغییرات مکانی لوژان حاصل از شبیه‌سازی گوسی به‌عنوان نمونه 
قو .سین آقرازهای. مات از ای هن شن (00 ضایی 


داده‌شده است. 


1) 


شکل ٩‏ (الف) نتایج بلوک سه‌بعدی شاخص کیفیت سنگ حاصل از شبیه‌سازی گاوسی در شبکه ایجادشده. (ب) نتایج شبیه‌ساز گوسی درترازهای مختلف با 
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شکل ۱۰- نقشه پلان مبانگین حاصل از ۱۰۰ شبیه‌سازی گاوسی لوژان در چند تراز ارتفاعی مختلف 


نشریهة مهندسی عمران فردوسی 


سال سی و چهارم. شمارة چهار: ۱۶۰۰ 


رضا دهشیبی - رضا جهانشاهی - مرنضی مظفری - امین عصاری 


در راستای ارزیابی عدم قطعیت مقادیر شبیه‌سازی‌شده 
متغیرهاء شبیه‌سازی‌ها بر روی یکدیگر منطبق شده و احتمال 
بزرگ یا کوچک‌تر از یک آستانه محاسبه‌شده است. با توجه به 
اينکه تشخیص زون‌های تراوا با مقادیر شاخص کیفیت سنگ 
پایین. مقادیر بالای لوژان و سیمان‌ خوری مدنظر است؛ نقشه‌های 
احتمال کوچک‌تر بودن مقادیر شاخحص کیفیت سنگ از مقدار 
حد آستانه ۳۰ درصد بالاتر بودن لوژان و سیمان‌خوری از مقادیر 
حد آستانه به ترتیب ۲ و ۱۰۰ کیلوگرم در واحد طول تهیه‌شده 
است. از این نقشه‌ها می‌توان در تعیین عدم قطعیت در راستای 
تعیین نقاط مناسب جهت تزریق سیمان و جلوگیری از فرار آب 
انظاده نود جهت شتاسایی نهتر تقاط بهیر اشک تفه عفر را 
در ترازهای مختلف ارتفاعی ترسیم و با یکدیگر مقایسه کرد. در 
اشکال (۱ ۱۲و (۱۳)ستد قبواته از نقشه احقبال:مقاوین مه 
متغیر شاخحص کیفیت سنگ. لوژان و سیمان‌خوری کوچک‌تر یا 
بزرگ‌تر از یک حد آستانه ترازهای مختلف مشاهده می‌شود. با 
توجه به این نقشه‌ها مشاهده می‌شود در محدوده‌ی مخزن و کمی 
جناح راست سد احتمال مقادیر شاخص کیفیت سنگ کوچک‌تر 
از ۳۰ احتمال بالاتر از ۲ و ۱۰۰ برای متغیرهای لوژان و 
سیمان‌خوری وجود دارد. به‌منظور بررسی رابطه بین گسل‌ها؛ 
نقشه‌های حاصل از میانگین احتمال شاعص کیفیت سنگ کمتر 
از با تشه گدا‌های سره ی مظقه نطین دادوشته ایخ 
(شکل ۱۶). نتایج نشان‌دهنده این است که در بدنه مخزن و جناح 
راست سد در نزدیکی و در محل تلاقی گسل‌ها احتمال فرار آب 
در اثر وجود مناطق با تراوایی بالا افزایش‌یافته است. اگر نقشه‌ی 
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حاصل از احتمال مقادیر لوژان بالاتر از ۲. نقشه احتمال شاخص 
کیفیت سنگ پایین‌تر از ۲۰ و نقشه احتمال سیمان خوری جهت 
مقادیر بزرگ‌تر ۱۰۰ کیلوگرم بر متر با یکدیگر مقایسه شود 
می‌توان دریافت که در دو محل مرکز مخزن سد و شمال غربی 
آن انطباق خوبی بین مقادیر پایین شاخص کیفیت سنگ. مقادیر 
بالای سیمان خوری و لوژان وجود دارد. درحالی‌که در بعضی از 
نقاط ازجمله نواحی شمال شرق مخزن. روند تغییرات مقادیر 
احتمال شاخص کیفیت سنگ و لوژان با یکدیگر سازگاری کمی 
نشان می‌دهند. البته از آنجاییکه همراه با برداشت داده های لوژان, 
شاخص کیفیت سنگ و سیمانخوری, اطلاعات لیتولوژی و درزه 
برداری ثبت نشده است؛ نمیتوان بطور دقیق مشخص کرد که آیا 
لوژانهای بالا یا شاخص کیفیت سنگ پایین مرتبط با سازنده‌های 
آهکی (کارستی) با درزه و شکاف های واحدهای سنگی هستند. 

قابلذکر است نتایج شبیه‌سازی هنگامی قابل‌قبول‌اند که قادر 
به دوباره‌سازی پارامترهای اماری و زمین‌اماری داده‌های اکتشافی 
باشند. بدین منظور نمودار فراوانی نمای به‌دست‌آمده از نتایج 
شبیه‌سازی ونمودارهیستوگرام داده‌های واقعی داده (داده‌های 
گمانه‌های اکتشافی) در نظر گرفته شد (شکل ۱۵). ملاحظه 
می‌شود که میانگین حاصله از نتایج شبیه‌سازی با میانگین حاصله 
از تحلیل آماری داده‌های گمانه‌های اکتشافی تقریباً برابر است 
درنتیجه نتایج حاصل از شبیه‌سازی به‌خوبی قادر به تولید مجدد 


پارامترهای آماری داده بوده است و از این لحاظ دارای اعتبار 


تنل 


شکل ۱۱- نقشه پلان احتمال متغیر شاخحص کیفیت سنگ با مقادیر کمتر از ۳۰ درصد (دوایر رنگی جهت نشان دادن مکان هایی که درآن ها متغییر مورد نظر 
مقادیر احتمال بالاتری را دارد) 


سال سی و چهارم» شمارة چهار: ۳۹ 
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شکل ۱۳ نقشه پلان احتمال سیمان‌خوری برای مقادیر با حد آستانه بالاتر از ۱۰۰کیلوگرم بر متر (دوایر رنگی جهت نشان دادن مکان هایی که درآن ها متغییر 
مورد نظر مقادیر احتمال بالاتری را دارد) 


شکل ۱۶ (الف) نمایش گسل‌ها در اطراف سد و (ب) نمایش مقادیر احتمال شاخحص کیفیت سنگ کمتر از ۳۰ بر روی سطح زمین 


نشریه مهندسی عمران فردوسی سال سی و چهارم. شمار؛ چهار» ۱۶۰۰ 


رضا دهشیبی - رضا جهانشاهی - مرنضی مظفری - امین عصاری 


شاخص کیفیت سنگ 
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شکل ۱۵ نمودار هیستوگرام مقادیر داده‌های اکتشافی برداشت از گمانه‌ها ومقادیرشبیه سازی شده 


در اين پژوهش مناطق دارای پتانسیل فرار آب درسدتنگاب با 
استفاده از شبیه‌سازی گاوسی متوالی موردبررسی قرار گرفت. با 
توجه به نقشه‌های شبیه‌سازی سه متغیر لوژان. شاخص کیفیت 
سنگ و سیمان‌خوری در دو محل واقع در مرکز مخزن و شمال 
خی ساسا وا ردنا تالا قتاما ی اس تقو ای 
مکان‌ها پتانسیل فرار آب وجود دارد. بررسی گسل‌های منطقه 
موردمطالعه نیز نشان داد تراکم گسل‌ها بیشتر در جناح راست سد 
بوده به همین دلیل وجود گسل‌ها می‌تواند سبب افزایش 
نفوذپذیری واحدهای سنگی شده باشد. بااین‌وجود در بعضی از 
نقاط به‌ویژه جناح چپ مخزن سد. ارتباط و همپوشانی بین مقادیر 
حداکثر لوژان و حداقل شاخص کیفیت سنگ. دیده نشده است 
که می‌تواند به دلیل ماهیت متفاوت و نحوه اندازه‌گیری دو متغیر 
شاخص کیفیت سنگ و وژان,باشد و همچنین یکسان نبودن عمق 
گمانه‌ها و برخورد گمانه‌ها در ترازهای مختلف ارتفاعی با 


مرخ 


لایه‌هایی با خصوصیات مختلف را نیز می‌توان از دلایل این 
موضوع ذکر کرد.درنهایت پيشنهاد می‌گردد تلفیق اطلاعات 
زمین آمار با اطلاعات ژئوالکتریک در شبیه‌سازی زمین آماری انجام 
شود و در مطالعات مشابه در هنگام اندازه گیری داده‌های شاخص 
کیفیت سنگ داده‌برداری درز و شکاف مغزه‌های سنگ انجام شود 
و مطالعات و تفسیرهای شاخص سنگ کیفیت بر اساس درز و 
شکاف صورت گیرد. همچنین باید اشاره کرد محل پیزومترها؛ 
ترریق ماده رنگی و برداشت نمونه های رنگی و هیدروشیمیایی 
در مطالعات قبلی در مناطق دورتر از موقعیت چاه های مشاهده 
ای و چاههای ژئوتکنیک کنونی قرار داشتند و در حال نیز آن 
پیزومترهای یا تخریب شده اند و یا سطح ایستابی آنها برداشت 
نمی شود. لذا در مطالعات پیشین» درمورد مسیر حرکت آب 
زیرزمینی و نقاط با نفوذپذیری بالا و یا کم در محدوده مطالعاتی 
ای تفگ بای خداشز انز که هران سقاییبه ابقر یوش 


۵۷۵۵ مناوتق۱ ۵۶ عم مامعونل متافقط‌هاو. و رمصمصمطتهملطمممظ ریت ر2صهاتهطاتااک ریظ رامومل بل بامنا20 .1 


2004(۰) ,751-7760 00۰ ,(1) 27 ۷۵۱۰ رکهع۵90۱0 ۷۷۵۵۲ 171 ۸0۷۵۱6۵5 


1969(۰) ,و۲۲۵ هت 0۵ ب۵۵ ار ومتاملاهاو0ع0) ۸۵۵11۵0 ربا یک به۷ه)9۳۱۷26 مباط .ظ م82 .2 


سال سی و چهارم. شمارة چهار» ۱۶۰۰ 


نشریه مهندسی عمران فردوسی 


۲۸ کاربرد زمی نآمار در شناسایی محل‌های با پتانسیل بالای فرا رآب در سد تنگاب 


0ص 10 میج صمع مصتااه0مصه وم فمتافتامادممع )صلممعامانهه ۵۴ ممتاه‌تانم ۳۸ ریخ رقعنام 1229521 رن رقصقصو تا بو 
(2009) ,1911 -1901 .00 ,(17)8 ۷۵۱۰ ,9۱۱۵ 60و۵0 ۰" (صیاتعآعظ) ععلننامه ۵60060-کو۵ه ۵ طا 01ح8صه 
15۵01 ۵۶ ممتاماصته. صهتوونادت لهتامم‌نومعو ها ناد ما ممتامعاع0 لناه۳. ریم رتتقطعفم بلط بتصقطحلو تصد۷۵ .4 
(2012) ,3737-3747 .00 ,(6)10 ۷۵۱۰ ,660562665 ۵ م0 ۸۳۵۱۵۲ ۲( ۴۲؟) ممتاهطع 1و6[ 
مصح تمعتاعتافاه عصتمد ملوعو له مط 2 عصتاملمصاهءنصطما۵ون. وبظ رتتهنط 1 رت ماما رثا رع8و۳07ظ بیط ,م۷272 .5 
2009(۰) ,67-76 00۰ ,(107)3-4 ۷۵۱ ,و000 ۵۱۱۵۵ ۳۲۵۵۵۵(۲۱) 6286 9926ع۳ عطل :واهها [2عتافتاهادهعع 
1 ساسانی نی ندا. بهرام حبیب نیاء و حسین قبادی. "تخمین تخلخل و توزیع فضایی تراوایی با استفاده از مدل زمین آماری دریکی از میدان‌های نفتی جنوب 
غربی ایران"» کنفرانس بین‌المللی علوم و مهندسی. امارت-دبی. موسسه ایده پرداز پایتخت ویرا؛ (۱۳۹۶). 
۷ عالی‌انوری» علی.» همایون کتیبه. و حمید محمودآبادی» "تخمین نفوذپذیری معادل در امتداد محور تونل امیرکبیر با استفاده از شبکه عصبی مصنوعی "" 
دومین کنفرانس مهندسی معدن تهران دانشگاه تهران (۱۳۸۷). 
۸ _ مشرفی فر» محمدرضا. کاظم برخورداری. و حمید موسوی. "ارزیابی خصوصیات ژئوتکنیکی و زمین‌شناسی مهندسی سد مخزنی ایوشان " دومین کنفرانس 
ملی مهندسی عمران. معماری و توسعه شهری بابل. مسسه علمی تحقیقاتی کومه علم آوران دانش؛ (۱۳۹۵). 
2 م1 «تاوتصمده بوصم ۵۶ ما2 مهد آقطمنعع۲ موه -کمتاملا0۵۵9)2. پنکا رمک منک ب7۷272 162901 وبا نان .9 
2016(۰) ,829 .00 ,75 ,۷۵۱ ,عهع۵6۱۵ ۵۲ ۳۱۱۷۱۳۵۱۱۱۹۱۵۵۸۵ ۱ عصتاعو 121660 ۳1 2762 ۲۵۵1 «تداصعمحصتل6و 
م10 منادیم ۵ ها فماه۷ 0۶عمنان 1 لصنه (نامک) ممتامصمنوع0 متلمتان 0۶1۵0۵1 وتوراقصش ری ر تمصصمطم]۷ ره رتتعفوفه .10 
۵ («60108) 0۳۱۵۱۱۵۵۲۵ ۱//۵1۲۷ظ و ص۱2 ,ماه جصدنا ۲۱۷ صیاتهک رطمههعممه صمت)هااصصته 141620 لمت)ممنا مهو عط معصتفت 
2016(۰) 1-12 .00 و(76)2 ,۷۵1 12۷۳۵۱۱۵۵ 1116 
متاترظ ه وه نصا و ممتاهام‌ووزنا ۷.9۵ رقاصمجصه‌اوعع] ۱۷ بز۴ ۶ ممتد‌ایصنه مصه ممتاهآهمعاما و0۵0 و۷۰ پذا ۱۷ 11۰ 
2010(۰) ,(20600۷۵۲ ۷) 12ماحصصیتامت 
۲ م لوا لمتاصمامم مملافمتصماومن متصموه 0۴ ممتاماله۳۲؟ رب ومهشنا متا بیط رعصع! مصتطاه بیط بتانا , ۷۷1982 .12 
(2004) ر173-188 00۰ و(321)1-3 ,۷۵1 ,۱۱۷۱۲۵۳۵۱۵۵ 01۵1 0۴1۵ 96۵6۵ (صه121۳) تعاننوه نامه ۷ 6ظ) م1 عصلعتن1۱0 
۳ عصاری, امین» و ضرغام محمدی. " شبیه سازی ضریب آبگذری باا ستفاده‌از روش‌های زمین‌آماری باهدف تعیین محل چاه‌های پمپاژ مطالعه موردی معدن 
گل‌گهر " اولین کنفرانس ملی مهندسی اکتشاف منابع زیرزمینی دانشگاه شاهرود. (۱۳۹۲). 
۲6۷۵ ۵6هام ۵ قفممتاهتنه ۷ 1۱۵8۲۵0۵۲۵ صرح آفمتاهمود ۵۶ فلوراممه آهتافتاماد۹060 ببظ رمتصصه‌طول0هی راو متفصطضه 14 
2007(۰) ,277-294 00۰ ,(1-3) 129 ,۷۵۱ ,۵و6و۸ ۵۱۱ ۱۷۵۱۱۱۱۵0۲۵ ۱۱۵۱۳۵0۵۵۵۵ 
۵ نخعی بهشته. " ستفاده از روش‌های شبیه سازی | ستوکا ستیک جهت تعیین زون‌های تراوا در سازند سخت پیت معدن گهر زمین. استان کرمان " ایران. 
بانال نامه کارشتاسی اوشت زمین شاسی باس حانشگام سیستان و بر ستان: ۱۳۹۷ 
۳ اسدی, ناصر» "بررسی مسئله فرار آب از سد تنگاب فیروزآباد باا ستفاده‌از ردیاب رنگی " پایان‌نامه کار شناسی ار شد آب شناسی-زمین شناسی, دانشگاه 
شیراز (۱۳۷۷). 
۶ 9006۵ معط ملهتامم00116 0 [۲۵0 2 28 دعصتنمه وتف ۵۶ ومتیمد مصصتا ۳۳/۹۱۵00۵16 رب ومتق2 و.ط مفلع2؟ رل متتطتتهک .17 
٩01192 ۷۵1۵۲ 2۵۵0۵15 ۵۵ ۲۵۵0115, ۷۵۱۰ 20)2(, 00. 138-147, )2005(۰‏ 
صمتماصته لمءتافتاماومعع ۵ عصتفناویی! 0۶ امعم ماع۷ع0 امتاهمه ۵۶ فلورلمصه پ«تفصتصصتاهزط م2 رتمصصم‌طم رب رتصمطلم .18 
2014(۰) ,1037-1047 00۰ ,(73)4 ,۷۵۱ 1۱۱۳۵0۱۱۵ ۱6 وه و6۵۵۵ 1۳۱۱۵۱۱۵۵۲۱۱۸۵ 0۴ ۱/1۱۷ 20۲020 
محمودی امین "محمد صدقی اصل. من صور پرویزی, و ار سلان مهرکی برر سی ذشت در جناح را ست سد خاکی مطالعه موردی سد خاکی تنگاب 


فیروزآباد " اولین کنفرانس ملی مهندسی ژئوتکنیک اردبیل دانشگاه محقق اردبیلی» (۱۳۹۲). 


نشریه مهندسی عمران فردوسی سال سی و چهارم. شمارة چهار» ۱۶۰۰ 


رضا دهشیبی - رضا جهانشاهی - مرتضی مطفری - امین عصاری 


۲۹ 


۰ کریمی حاجی.. عزت... رئیسی و محمد زارع آبزرشی, رابطه بین تغییرات سطح آب نفوذپذیری لابه‌های زمین‌شناسی بااستفاده‌از شبکه پیزومتری " اولین 


کنفرانس زمین‌شناسی مهندسی و محبط زیست تهران. انجمن زمین‌شناسی مهندسی ایران. دانشگاه تربیت‌معلم. (۱۳۷۸). 


۱ جاور <سن. و علی‌محمد ح<سین نژاده "آب‌بندی سد تنگاب فیروزاباد فارس با خروج مواد ریزدانه از شکیلات ر سوبی و سنگی و طراحی گمانه‌های 


متقاطع " دومین کنفرانس ملی نیروگاه‌های آبی کشور شرکت توسعه منابع آب و نیروی ایران» (۱۳۸۷). 

۲ 10۲ ۸۵061 13067 ۲6[ قفهرنا ۷۱۵011160 رت ب0عتفداظ وی بت۷]۵0 یط بل۷]۵0 ب.ه بتتامهتط2 ۷ب رصفلهکبا0ظ 
۵ , (۲12عع۸1) ههد صتفام ۱۷۵0۵ مطا صا قر) مه 08مطاعصه متافتافاوهعع مصئوت عصاممممه. وانازماهتمصاناه 
۰ ۰ ,(2021) ,۳0۲۵755 
:۰ ۷/۱ موجه لهمتامواهاوممم ماملنه ۷۲۵۱0 رخ ۷۲۰ رتصحفعظ تت2ع۸ .۲ .۲ هدومن ماما بهعطاصطژهن وبلط .ظ بهات۲2 
2021(۰) .174-184 00 ,(130)3 ۷۵۱۰ ,5661۱66 ۲۵۳۸ ۸۳۵/۱20 ر اصمحصم‌ودعهه اصتهاععصه 10۲ ممتاه‌نام مه هه 
6 1 طمتاه1۳)6۵01 المکصته: ۲۵۲ و0مطاممهر متافتمتماعل مه اهتافتاماو0۵0 ربع روم‌ه۷]۵0 ر.ظ تصهتزظ رتالاهژ 
2020(۰) .2678-2692 و0 ب(65)16 ,۷۵1 ,]02۲۱۵ 566665 ۳۵۲۵/۵96۵ و ص1۳2 رصتفد ۱۱0 
۵ 2 0۷ 0رهه متصهوه لزمو ۵ ممتامتلمين مه ممتممه ,رل ملتتهط .تا مهتطفن۷ و. مکلتتط«مطظ بو بکلناله]1۷ 
2020(۰) . ملقطمتاحصعاص1 .. ماتجهمعن .. ر هلق صتحعها .. لصه .. مصتفمعه .. مامصصع: .. هصلفدا.. متمتصطهعا.. لمعتافتاماد0عع 
1 106 .0 001.012/1 
1 01 ولوراممح مه عممممهه آه۲ممجم-9۵200 ریک .نا رقطتقطا ریکلیگ بتک م۱ بقط۵ع]۱۷ ریگ رتقصناک ر.ظ ,ت22 
0 2۵4 ماهنه 0اه رآ(ظ) ممتامن4صا متامصمممهما0ع۱ظ طمبامعطا طمهمممه ‏ تعنص مه :اتطنلهه 
106049.0۰ 10۰1080/10 /101.012 .(2021) رلقطمتامصمام ماتهعمم رد ومناوتصطهما 
معا متام ۲۵۵۲۵8۵0)20۷۵ ۵۶ طمتامصتصهامنا ری بتاکر ررظ مصقطا و2 م۷1 بت) م20 آ و ,9018 ب.9 ,۷۷228 ,۷۷.۰ ,21202 
از ۵00 01/۱۷۳۵۵۱۱۵ 0۵۵ ۱/0۵6۵ طمه200۲0 ومتافتافادههع مصلوتا 0]۲00۱فنصج تحعصنا ۵ ممامع0تفصمن )۷1 
189.2020۰ 080/19648 101.012/10.1 2)0200(۰) ,۳۲۱8۵1166۳۱۱۱8 
1989(۰) ,وفع وه ۱۲ ۵0 0و۱ روم‌تافتلهادهع0 ۸۵۵1160 رب یک به۷ه)29 9۲۱۷ وبلط .ظ ,1552 

.)۱۳۸۶( سسسنی پاک» علی‌اصغر» محمد شرف‌الدین» "تحلیل داده‌های اکتشافی " انتشارات دانشگاه تهران‎ ٩ 


۰ مدنی. حسن. "مبانی زمین‌آمار ‏ دانشگاه صنعتی امیرکبیر تهران؛ (۱۳۹۱). 
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سال سی و چهارم. شمارة چهار» ۱۶۰۰ نشريهة مهندسی عمران فردوسی 


۳۰ کاربرد زمی نآمار در شناسایی محل‌های با پتانسیل بالای فرا رآب در سد تنکاب 


نشریه مهندسی عمران فردوسی سال سی و چهارم. شمارة چهار» ۱۶۰۰ 


